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ABSTRACT 

, We calculate axisymmetric oscillations of rotating neutron stars composed of the sur- 

■ face fluid ocean, solid crust, and fluid core, taking account of a dipole magnetic field 

as strong as B$ ~ 10 15 G at the surface. The adiabatic oscillation equations for the 
solid crust threaded by a dipole magnetic field are derived in Newtonian dynamics, 

I , on the assumption that the axis of rotation is aligned with the magnetic axis so that 

■ perturbations on the equilibrium can be represented by series expansions in terms of 
CO \ spherical harmonic functions VJ m (r?, <f>) with different degrees / for a given azimuthal 

wave number m around the the magnetic axis. Although the three component models 
| can support a rich variety of oscillation modes, axisymmetric (to = 0) toroidal jt„ and 

■ spheroidal is n shear waves propagating in the solid crust are our main concerns, where 
\Q \ I and n denote the harmonic degree and the radial order of the modes, respectively. In 

■ the absence of rotation, axisymmetric spheroidal and toroidal modes are completely 
\ decoupled, and we consider the effects of rotation on the oscillation modes only in 

O ^ the limit of slow rotation. We find that the oscillation frequencies of the fundamental 

^ \ toroidal torsional modes it n in the crust are hardly affected by the magnetic field as 

j_i ■ strong as B$ ~ 10 15 G at the surface. As the radial order n of the shear modes in the 

\ crust becomes higher, however, both spheroidal and toroidal modes become suscepti- 

C$ ■ ble to the magnetic field and their frequencies in general get higher with increasing 

\ Bs- We also find that the surface g modes and the crust/ocean interfacial modes are 

V~j ■ suppressed by a strong magnetic field, and that there appear magnetic modes in the 

r> \ presence of a strong magnetic field. 
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1 INTRODUCTION 

Recent discoveries of quasi-periodic oscillations (QPOs) in the giant flares of Soft Gamma-Ray Repeaters SGR 1806-20 (Israel 
et al 2005) and SGR 1900+14 (Strohmayer & Watts 2005) and the confirmation of the discoveries (Watts & Strohmayer 2006, 
Strohmayer & Watts 2006) have made promising asteroseismology for magnetors, neutron stars with an extremely strong 
magnetic field (see, e.g., Woods & Thompson 2004 for a review on SGRs). Israel et al (2005) discovered QPOs of frequencies 
~ 18, ~ 30 and ~90Hz in the tail of the SGR 1806-20 hyperflare observed December 2004, and suggested that the 30Hz 
and 90Hz QPOs could be caused by seismic vibrations of the neutron star crust (see, e.g., Duncan 1998). Later on, in the 
hyperflare of SGR 1900+14, Strohmayer & Watts (2005) found QPOs of frequencies 28, 53.5, 84, and 155 Hz, and claimed that 
the QPOs could be identified as low I fundamental toroidal torsional modes in the solid crust of the neutron star. Although 
the interpretation in terms of the crustal torsional modes is promising, the mode identification cannot always be definite. 
In fact, if different classes of oscillation modes can generate similar periodicities and information other than the periods are 
not available for the QPOs, it is usually difficult to assign a class of oscillation modes to the periods observed in preference 
to other classes of modes. This is particularly true for high frequency QPOs (e.g., 625Hz QPO, Watts & Strohmayer 2006; 
1835Hz QPO and less significant QPOs at 720 and 2384 Hz in SGR 1806-20, Strohmayer & Watts 2006), for which there exist 
various classes of modes that can generate the periodicities observed. In this case, the pattern of observed frequencies, that 
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is, the observed frequency spectrum could be a key for mode identification and hence understanding the underlying neutron 
stars themselves. McDermott et al (1988) have carried out detailed modal analyses for neutron star models with a solid crust, 
but without the effects of magnetic field and rotation on the oscillation modes. Since the detected QPOs are thought to come 
from magnetors inferred to possess a strong magnetic field, it is useful to calculate oscillations of rotating neutron stars that 
have a solid crust and a strong magnetic field. 

An extensive theoretical modal analysis of magnetic neutron stars having a solid crust was first carried out by Carroll et 
al (1986), who however employed a cylindrical geometry for the analysis of the neutron star models, and assumed a uniform 
magnetic field whose axis is parallel to the axis of the cylinder. The maximum strength of the magnetic field they examined is 
Bs ~ 10 12 G, which may be too weak for mangetors for which Bs >10 14 G is inferred. Carroll et al (1986) have suggested the 
existence of Alfven modes and also found the transformation of the g modes in the surface ocean into magnetic modes as the 
field becomes stronger. More recently, Piro (2005) solved a simplified set of oscillation equations for toroidal torsional modes 
employing improved shear modulus (Strohmayer et al 1991) and up to date microphysics for equations of state to construct 
back ground neutron star models. Since the fluid core was ignored in the calculations by Carroll et al (1986) and Piro (2005), 
no reliable modal analyses were possible for spheroidal modes that can have substantial amplitudes in the core. In this paper, 
we calculate various oscillation modes of neutron star models that have a solid crust and are threaded by a dipole magnetic 
field. Since a non-radial mode of a spherical neutron star thereaded by a dipole magnetic field cannot be represented by a 
single spherical harmonic function, we employ series expansions to represent the perturbations accompanied by a mode. In 
§2, we give a brief description of method of solution we empoy, numerical results are given in §3, and §4 is for conclusions. 
The oscillation equations are given in Appendix as well as boundary and jump conditions used in this paper. 



2 METHOD OF SOLUTION 
2.1 Perturbation Equations 

To derive oscillation equations for a magnetized and rotating neutron star with a solid crust, we follow Carroll et al (1986), 
Lee & Strohmayer (1996), and Lee (2004). We consider no general relativistic effects on oscillations, that is, the oscillation 
equations are derived in Newtonian dynamics. We employ spherical polar coordinates (r,9,<f>), whose origin is at the center 
of the star and the axis of rotation is the axis of 9 — 0. We assume a dipole magnetic field given by 

B = /i m V(cos #/r 2 ), (1) 

where p m is the magnetic dipole moment. For simplicity, we also assume that the magnetic axis coincides with the rotation 
axis. Since the dipole field is a force-free field such that (V x B) x B = 0, the field does not influence the equilibrium structure 
of the star. Assuming the axis of rotation is also the magnetic axis, the temporal and angular dependence of perturbations 
can be given by a single factor e i < m * +a '*> , where m is the azimuthal wavenumber around the rotation axis and to = a + mQ is 
the oscillation frequency observed in the corotating frame of the star where a is the oscillations frequency in an inertial frame, 
and fl is the angular frequency of rotation. The linearized basic equations applied in the solid crustal region of the star are 
then given by 

- u; 2 £ + 2iu>il x £ = -V ■ a' - ■ a + — !— (V x B') x B -^—iujE' x B, (2) 
P + V • {pQ = 0, (3) 

B' = V x (£ x B), (5) 

E' = x B/c, (6) 

where p is the mass density, p is the pressure, c is the velocity of light, £ is the displacement vector, B' and E' are the Euler 
perturbations of magnetic and electric fields, respectively, and the other physical quantities with a prime (') denote their Euler 
perturbations, and A is the Schwartzshild discriminant defined by 

_ rilnp 1 dlnp 
and 

r-fel • (8) 



dlnp 
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Note that we have employed the Cowling approximation neglecting the Eulerian perturbation of the gravitational potential, 
and that no effects of rotational deformation are included. In equation (2), er' denotes the Euler perturbation of the stress 
tensor and is obtained from the Lagrangian perturbation defined in Cartesian coordinates by 

&(Jij = (F 1 pu)S i j + 2/u(wij - ^uSij) (9) 
with Uij being the strain tensor defined by 

i (db , di 3 



" + < 10 > 

where Sij denotes Kronecker delta, /i is the shear modulus, and u = Xw-i Un - The last term on the right-hand side of equation 
(2) represents the contribution from the displacement current, where infinite conductivity has been assumed. The linearized 
basic equations for a fluid region may be obtained by simply replacing the terms V • a and V • <r' by — Vp and — Vp', 
respectively. 

Since the angular dependence of perturbations on a rotating and magnetized star cannot be represented by a single 
spherical harmonic function, we expand the perturbed quantities in terms of spherical harmonic functions Y t m with different 
Is for a given m, on the assumption that the axis of rotation coincides with that of the magnetic field. The displacement 
vector £ and the perturbed magnetic field B' are then approximately represented by finite series expansions of length j max as 



3 = 1 

and 
B' 



{[S h (r) + H h (r)V] Y%(6, 0) + T t , (r) e r x VY^O, 0)} e™\ (If) 



Jmax 

£ {[&jj(r) + &/?(r)v] Y l T{e,4>) + bJ j {T) e r xVl^(M)}e w , (12) 



B (r) 

and the pressure perturbation, p', for example, is given by 



p' = 2pi i (r)l£ , (M)e iu ' t . (13) 
j'=i 

where 50(5*) = fJ-m/r 3 , and ij = |m| + 2(j — I) and lj = lj + 1 for even modes, and = |m| + 2 j — f and = lj — 1 for odd 
modes, respectively, and j = 1, 2, 3, • ■ ■ , j max . In this paper, we have used jWx = 12. Substituting the expansions (11)~(13) 
into the linearized basic equations (2)~(6), and making use of equations (9) and (10), we obtain a finite set of coupled linear 
ordinary differential equations for the expansion coefficients such as SV (r) and bf, (r), which we call the oscillation equations 

3 

solved in the solid crust. The oscillation equations to be solved in magnetic fluid regions can be derived in the same manner 
(e.g., Lee 2004). The sets of oscillation equations thus obtained for the solid and fluid regions are given in Appendix A for 
the case of axisymmetric modes with m = 0. The boundary conditions at the center and the surface of the star and the jump 
conditions imposed at the interfaces between fluid and crustal regions are discussed in Appendix B. It is important to note 
that in the case of Q = the oscillation equations for axisymmetric modes with m = are decoupled into those for spheroidal 
modes and toroidal modes, respectively, and that axisymmetric spheroidal and toroidal modes are coupled only through the 
effects of rotation. For non-axisymmetric modes with m =fc 0, however there occurs no decoupling between spheroidal modes 
and toroidal modes even for £7 = 0. 



2.2 Non-Magnetic Core 

In a fluid region where the magnetic pressure p_g = B ( )(r)/87r is much smaller than the gas pressure p, that is, where the 
Alfven velocity va = ^2ps/p is much smaller than the sound velocity ~ \Jp/ p or rui, the magnetic perturbations b H and 
ib T suffer very rapid spatial oscillations except for extremely low frequency oscillations (e.g., Lee 2004, see also Biront et al 
f982, Roberts & Soward 1983, Cambell & Papaloizou 1986). This rapid spatial oscillation happens in the fluid core for both 
spheroidal and toroidal modes. For axisymmetric toroidal modes, for example, the dimensionless wavenumber k r in the radial 
direction may be given by 

k r ~ ruj/vA, (14) 

and the wavenumber can be very large in the magnetic fluid region immediately below the crust where ps/p <S 1 even for 
Bs = Bo(R) ~ I0 15 G at the surface. Since this kind of rapid spatial oscillations usually do not occur in the solid crust, 
there exists an abrupt change in the property of the magnetic eigenfunctions across the boundary between the solid crust 
and the fluid core. Numerically, this means that we need an extremely large number of mesh points to correctly calculate 
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very short magnetic perturbations in the fluid core, which is not feasible. Since the rapid spatial oscillations in the magnetic 
perturbations come from the fact that the magnetic pressure significantly weaker compared to the gas pressure in the fluid 
core, it is legitimate to assume that the fluid core is non-magnetic so that the oscillation equations for non-magnetic fluid 
stars (e.g., Lee & Saio 1990) could be used. Assuming non-magnetic core, we neglect a possible dissipation of the oscillation 
energy due to Joule heating in the core. 



3 NUMERICAL RESULTS 



In this paper we mainly discuss the case of ft = 0, in which axisymmetric (m = 0) toroidal and spheroidal modes are decoupled. 
We briefly discuss the effects of rotation on the axisymmetric modes, assuming slow rotation so that 
which case toroidal and spheroidal modes are only weakly coupled. 



n/^/GM/R s 



« 1, in 



3.1 Equilibrium Models 

The neutron star models we use consist of the surface fluid ocean, the solid crust, and the fluid core and are the same as those 
employed by McDermott et al (1988) for modal calculations without rotation and magnetic fields. The models are obtained 
from the fully general relativistic evolutionary cooling calculations of neutron stars by Richardson et al (1980). The outer 
crust extends down to the neutron drip point at p = 4.3 x 10 1 g cm and is assumed to consist of bare Fe nuclei embedded 
in a uniform, neutralizing, degenerate electron gas. The strength of the Coulomb interaction between ions is characterized by 
the dimensionless parameter V, the ratio of the Coulomb energy to the termal energy ksT with ks being Boltzmann constant. 
The matter of the outer crust is assumed to undergo a first order fluid/solid phase transition at F = 155, and there exists 
the fluid ocean above the crystallization boundary at F = 155. In addition to the thermodynamic contributions from the 
nuclear and electronic kinetic energies, the equation of state used in the outer crust includes contributions from photons, the 
Coulomb interactions between nuclei and electrons, and nuclear vibrations and rotations. The inner crust extends from the 
neutron drip point at 4.3 x 10 11 g cm~ 3 to the base of the crust at 2.4 x 10 14 g cm~ 3 and it is assumed to consist of nuclei 
with Z ~ 40, degenerate electrons, and degenerate, nonrelativistic neutrons, where the actual composition of the nuclei in 
the inner crust is taken from Negele & Vautherin (1973). For the inner crust, the zero temperature equation of state is that 
by Negele & Vautherin (1973), and the leading order thermal corrections for the nuclei, electrons, and free neutrons are also 
included. At densities greater than 2.4 x 10 14 g cm" 3 , the lattice is assumed to dissolve, and the core of the neutron star is 
taken to consist of a mixture of free and highly degenerate neutrons, protons, and electrons. The equation of state in the fluid 
core is that by Baym, Bethe, & Pethick (1971), and are added the leading thermal corrections for the three species. The outer 
most fluid envelope models are from the calculations by Gudmundsson et al (1983), and are fitted to the evolutionaly core 
models by Richrdson et al (1980) at mass density p w 10 10 g cm -3 . The more details of the models for oscillation calculations 
can be found in McDermott, van Horn, & Hansen (1988). 

The shear modulus p of the solid lattice is that given by Pandharipande, Pines, & Smith (1976): 

/* = 0.3711 21/3 N , (15) 

where tin is the number density of the nuclei. An improved calculation for the shear modulus was given by Strohmayer et 
al (1991). This imporved shear modulus was used by several authors to calculate crustal oscillations of neutron stars, for 
example, by Duncun (1998) and Piro (2005). In this paper, we did not attempt to use imporved models for modal calculations 
since our main concerns are to investigate the modal properties of magnetic neutron stars and not to try a detailed comparison 
between theoretical predictions and observational results for the QPOs from the SGRs. 

For mode calculations in this paper, we mainly use neutron star models named NS05T7 and NS13T8, and the mass 
and radius of the former are M — O.5O3M0 and R = 9.839km, respectively, and those of the latter are M — 1.326Mq and 
R — 7.853km, respectively. Both of the models have a solid crust, and the thickness of the crust is Ar/R ~ 0.24 for NS05T7 
and Ar/R ~ 0.055 for NS13T8. The thickness of the surface ocean, on the other hand, is Ar/R ~ 3.8 x 10~ 5 for NS05T7 and 
Ar/R ~ 2.3 x 10~ 3 for NS13T8, respectively. The details of the physical parameters of these models are given in McDermott 
et al (1988). 



3.2 Toroidal Modes 

In Figure 1, the normalized frequencies Q = u/ 'v/ GM/R 3 of the axisymmetric (m = 0) toroidal torsional modes vt n are 
plotted versus Bs for the models NS05T7 (top panel) and NS13T8 (bottom panel), where Bs = p-m/R is the strength 
of the dipole magnetic field at the surface. The frequencies Co are local ones and for convenience the redshifted frequencies 
defined by v x = to(l — 2GM / Rc 2 ) 1 ^ 2 /2n are also displayed on the right axis of the panels. We find that the frequencies Co 
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of the fundamental toroidal modes are only weakly affected by the dipole magnetic field of strength as large as Bs ~ 10 15 G 
at the surface, which is consistent with the results by Piro (2005). Note that the frequencies of the fundamental modes yto 
show a slight decrease with increasing Bs for low values of I' . The overtone torsional modes yt n having n > 1 become more 
susceptible to the magnetic field as the radial order n gets higher and the wavelengths in the radial direction become shorter 
so that several wavelengths of the modes are spaned by the strong magnetic region in the outer most envelope (or crust). We 
find as a general trend the frequency of the overtone modes increases with increasing Bs- Besides the magnetic effect just 
mentioned, there exists a different kind of magnetic effects closely related to the fact that the separation of variables using 
a single spherical harmonic function is not possible for the oscillations of stars with a dipole magnetic field. In this case we 
may consider crustal toroidal modes with different I's but with the same parity are coupled in the presence of the magnetic 
field. Since the toroidal modes with a given radial order n > 1 are nearly degenerate in frequency for different values of the 
coupling effect brings about interference between them, which becomes significant as the field strength increases. As discussed 
by McDermott et al (1988) for non-magnetized neutron stars, the frequencies of the fundamental toroidal modes with different 
I's approximately scale as 1 



which is also confirmed even in the presence of a strong magnetic field Bs ~ 10 G. 

The eigenfunctions iTy and ibf of the axisymmetric fundamental 2to mode of the model NS05T7 are depicted in Figure 
2 for Bs = 10 12 G in panels (a) and (b) and for Bs = 10 15 G in panels (c) and (d), where the solid and dashed lines stand for 
/' = 2 and 4 for iTy and for 1 = 1 and 3 for ibf ', and the amplitude normalization is given by iTy =2 = 1 at the surface of the 
star. Different from the case with no magnetic field, the toroidal components iTy of the displacement vector is continuous at 
the crust/ocean interface (see Carroll et al 1986). However, since we have assumed the fluid core is non-magnetic, no toroidal 
components iTy and ibf exist in the core for Q = 0. Both for Bs = 10 12 G and Bs = 10 15 G, iTy =2 dominates iTy =A in 
the inner crust where pb <C p. For Bs = 10 12 G, the eigenfunctions iTy and ibf are approximately constant as functions of 
r/R, but for Bs = 10 15 G they are influenced by the strong magnetic field in the outer crust where the magnetic pressure pb 
dominates the gas pressure p. 

Reflecting the difference in the thickness of the solid crust between the two neutron star models, the normalized frequency 
spectra of the toroidal torsional modes show substantial differences. For example, the fundamental yto modes of the model 
NS13T8 have normalized frequencies u> lower than those of the model NS05T7, while the overtones of the former have higher uj 
than the latter. This property remains the same even in the presence of a strong magnetic field. For the periods of the torsional 
modes of non-magnetic neutron stars, McDermott et al (1988) gave extrapolation formulae, in which the local periods are 
proportional to the radius of the star for the fundamental modes and to the crust thickness for the overtone modes for a given 
spherical harmonic degree I' . The normalized frequencies u> in the figure are consistent with what the extrapolation formulae 
predict. 

On a closer look at the behavior of the low frequency torsional modes of the model NS13T8 in Figure 1, one may find 
that the torsional modes suffer mode crossings with magnetic modes, for which the oscillation frequency increases rapidly 
with increasing Bs and the oscillation energy is dominated by the magnetic perturbations. In fact, if we calculate the energies 
defined by 



we find that Eb(R) is much larger than Ek{R) for the magnetic modes. Note that V • £ = for toroidal magnetic modes. For 
more detail, see below §3.3. 

3.3 Spheroidal Modes 

Figure 3 plots the oscillation frequencies lj of the low radial order spheroidal shear modes ; = 2S n and the core/crust interfacial 
modes ii 2 for the models NS05T7 (top panel) and NS13T8 (bottom panel). For the model NS05T7, also plotted are magnetic 
modes, labeled rrik, that are found in the presence of a strong magnetic field. Note that in this paper, we use the symbols iii and 
ii2 to denote the crust/ocean and core/crust interfacial modes, respectively, where the amplitudes of the former are strongly 
localized in a narrow region around the crust /ocean interface and the amplitudes of the latter are largest at the core/crust 
interface. The oscillation frequencies w of the spheroidal shear modes increase with increasing Bs and this frequency increase 
becomes more rapid for higher radial order n modes. As in the case of toroidal torsional modes, for a given radial order n, 
the spheroidal shear modes is n with different degrees I are nearly degenerate in frequency, and mode couplings between the 

1 As pointed out by one of the referees, the frequency scaling formula given by equation (16) is correct only for /' 2> 1, and a much 



better scaling formula should be u>y ~ \J (I' — + 2)/4 uiy =2 as one can easily verify by solving the simple case of a uniform density, 
uniform shear modulus oscillating crust. 



y/V{V + l)/6 u v 



=2, 



(16) 




(17) 
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shear modes become significant for strong Bs- Figure 3 also shows that the frequency of the core/crust interfacial modes 112 
is hardly modified by the magnetic field as strong as Bs = 10 15 G, except that the interfacial modes of the model NS05T7 
expreience avoided crossings with magnetic modes, whose oscillation frequencies rapidly increase with increasing Bs- 

As in the case of the toroidal torsional modes, the normalized frequency spectra of the spheroidal modes of the two models 
appear differently, reflecting the difference in the thickness of the crust. The different appearance for the shear modes of the 
two models can be understood by using an approximation formula for the period, which is proportional to the crust thickness 
(e.g., McDermott et al 1988). 

The eigenfunctions Hi and bfi of a 2S1 mode of the model NS05T7 are plotted in Figure 4 for Bs = 10 12 G in panels (a) 
and (b) and for Bs = 10 15 G in panels (c) and (d), where the amplitude normalization is given by Si t (R) = 1 at the surface. 
Note that because of the jump conditions applied at the interfaces between the solid crust and the fluid zones (see Appendix 
B), the horizontal component of the displacement vector is continuous at the outer interface, but it is discontinuous at the 
inner interface, at which almost free slippery jump condtions are employed. For Bs = 10 12 G, the eigenfunctions bf? exhibit 
zigzag behavior in the outer crust reflecting that of Ti, which is caused by discontinuous change of equilibrium composition 
with increasing p. For Bs = 10 12 G, the expansion coefficients H tl and bp are dominant over others, but for Bs = 10 15 G the 

first few components of Hi and bff are comparable to each other, and the eigenfunctions in the outer crust are significantly 
modified by the strong field comaperd with those for the case of B s = 10 12 G. As shown by Fi gure 4, the amplitude of b v 
for Bs = 10 12 G is much larger than that for Bs = 10 15 G for the same amplitude normalization Si 1 (R) = 1. However, if 
we compare the quantities Ek(t) and i?s(r), we find Ek(R) is much larger than Eb(R) for Bs = 10 12 G, but Ek(R) and 
Eb(R) are comparable to each other for Bs = 10 G so that the magnetic perturbations is important to determine the modal 
properties of the shear modes in the presence of a strong magnetic field. 

In Figure 5, the frequency Co of the fundamental 2/ mode of the model NS05T7 is plotted versus Bs- The 2/ mode heavily 
suffers mode crossings with high radial order shear modes is n as Bs increases, and no pure 2/ mode could be identified in the 
presence of a strong magnetic field. This is also the case for high frequency p modes for a strongly magnetized star. Note that, 
even at Bs ~ 10 12 G, the fundamental mode 2/ with uj = 1.8862 is affected by the shear mode 2S n =io that has the frequency 
Q = 1.8726 at B s = 0. 

The frequency u of the <?i mode in the fluid ocean and that of the core/ocean interfacial mode ; = 2«i are plotted versus 
Bs for the model NS05T7 in Figure 6. The frequencies begin to decrease rapidly to zero beyond Bs ~ 10 5 G, suggesting the 
strong magnetic field suppresses these modes. This suppression may come from the dominance of the magnetic force over 
the buoyant force in the fluid ocean at large Bs- We obtain almost the same suppression of the surface g modes and the 
crust/ocean interfacial modes for the model MS13T8, but they can survive much stronger mag netic field B s ~ 10 10 G. This 
may be because the fluid ocean of the hotter model NS13T8 is much thicker and has a more developed radiative region, 
compared to the model NS05T7, giving much stronger buoyant force for the modes. It is interesting to note that the result 
for the g and iii modes in the fluid ocean in this paper is different from that obtained by Carroll et al (1986), who suggested 
that the frequency of gi proportionally increases with increasing Bs to be a magnetically dominating mode, which they called 
g/m modes. The reason for the discrepancy is probably the difference in the geometry of the models and the magnetic field 
configulation between the two calculations. 



3.4 Magnetic Modes 

For both toroidal and spheroidal modes, we find oscillation modes that can be regarded as magnetic modes, for which the 
frequency Co rapidly increases with increasing Bs, and the magnetic energy Eb(R) is dominant over the kinetic one Ek(R), 
although its appearance looks different depending on the model structure, particularly, on the crust thickness. It is also 
important to note that there seem to exist critical strengths of the magnetic field beyong which magnetic modes are allowed 
to exist. We find that the low frequency torsional modes of the model NS13T8 are often affected by mode crossing with the 
magnetic modes, but that those of the model NS05T7 rarely suffer from the mode crossings within the frequency and magnetic 
field strength ranges we investigated in this paper. In the same parameter ranges as those for the toroidal torsional modes, 
we find magnetic modes interacting with the core/crust interfacial modes ;i2 for the model NS05T7, but no examples of the 
mode crossings of this kind are found for the model NS13T8. 

As an example of the mode crossings between the torsional modes and magnetic modes found for the model NS13T8, 
we plot, in Figure 7, Co of the 2^0 (left panel) and zti (right panel) modes versus log_B s . The panels are magnifications of 
the corresponding parts from Figure 1, and clearly indicate mode interactions with magnetic modes as Bs increases. It is 
interesting to note that between Bs = 10 12 G and ~ 10 13 G, the 2*1 mode also interact with sequences of modes whose 
frequencies decrease as Bs increases. Figure 8 shows an example of the eigenfunctions of a toroidal magnetic mode having 
Co — 0.01557 at Bs = 10 14,48 G for NS13T8, where the amplitude normalization is given by iTy = 2 = 1 at the surface. For this 
mode, the expansion coefficients iT;/ =2 and ibj =1 are dominant, and the magnetic perturbation ibj shows almost discontinuous 
change at the crust /ocean interface, and has almost negligible amplitudes in the fluid ocean. Figure 9 is an example of the 
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Table 1. Coefficient C 2 for axisymmetric (m = 0) modes for the models NS05T7 and NS13T8 







IN OUO -L t 






1\TG1 qTS 




Mode 




Period (ms) 


c 2 




Period (ms) 


c 2 


2*0 


0.04019 


18.6 


-6.75 x 10- 2 


0.01897 


17.4 


-1.26 x lO" 1 


2*1 


0.3286 


2.29 


2.76 x 10 1 


0.4534 


0.727 


2.10 x 10 2 


3*0 


0.06352 


11.8 


-6.20 x 10~ 2 


0.03000 


11.0 


-8.31 x 10~ 2 


3*1 


0.3318 


2.26 


1.31 x 10 1 


0.4539 


0.726 


1.05 x 10 2 


4*0 


0.08517 


8.81 


-9.83 x 10~ 2 


0.04024 


8.19 


-6.39 x 10~ 2 


4*1 


0.3361 


2.23 


6.74 x 10° 


0.4546 


0.725 


6.94 x 10 1 


2H 


0.008282 


90.6 


1.03 x 10 2 


0.03326 


9.91 


2.56 x 10 1 


3H 


0.01171 


64.1 


7.95 x 10 1 


0.04703 


7.01 


1.92 x 10 1 


2*2 


0.1029 


7.30 


1.10 x 10 1 


0.01862 


17.7 


1.51 x 10° 


3*2 


0.1400 


5.36 


8.33 x 10° 


0.02889 


11.4 


2.52 x 10° 


2 Sl 


0.3093 


2.43 


-2.10 x 10 1 


0.4512 


0.731 


-1.92 x 10 2 


2S2 


0.5556 


1.35 


-5.37 x 10 1 


0.7704 


0.428 


-1.84 x 10 2 


3S1 


0.2884 


2.60 


-2.10 x 10 1 


0.4494 


0.734 


-1.03 x 10 2 


3S2 


0.5518 


1.36 


-9.51 x 10° 


0.7684 


0.429 


-1.02 x 10 2 


2/ 


1.886 


0.398 


3.11 x lO" 1 


1.434 


0.230 


1.39 x 10-1 


2P1 


4.038 


0.186 


2.09 x 10- 1 


3.967 


0.0831 


1.42 x 10~ 2 



eigenfunctions of a spheroidal magnetic mode of Co = 0.07965 at Bs = 10 1 ' G for the model NS05T7, where the amplitude 
normalization is given by Si=o(R) = 1. For this magnetic mode the first components Hj_ 2 and bff =1 of the series expansion 
are not necessarily dominant over other components associated with higher Is. Note that the relative thickness of the surface 
ocean of the model NS13T8 is much larger than that of NS05T7. 

3.5 Effects of Slow Rotation 

For slow rotation, we may approximate the oscillation frequency of axisymmetric (m = 0) modes as 

Ld^LJO + GW 2 , (18) 

where u>o denotes the oscillation frequency for the non-magnetic and non-rotating case. Note that the linear term in 
is proportional to m and does not appear for axisymmetric modes, and that no effects of rotational deformation of the 
equilibrium structure on the oscillation are included in our calculation and hence the second order effects exclusively come 
from the Coriolis force, which is dominant for low frequency modes with Q <1 (e.g., Lee 1993). We tabulate in Table 1 the 
coefficients C2 for several toroidal and spheroidal modes including torsional ;/*„ and shear is n modes for the models NS05T7 
and NS13T8. Note that we have used the labels iii and 112 to denote respectively the crust/ocean and core/crust interfacial 
modes for both the models. For the model NS05T7, the frequency of the 2*1 mode for Bs — is lower than that of the 
212, but for the model NS13T8 the frequency of the 2*1 mode is higher than that of the 212 mode. From the table, we find 
that the second order response of the fundamental toroidal modes to rotation is negative and small, but that of the \it n =\ is 
positive and rather large in the sense that the frequency correction C2& 2 can be comparable to (Do itself for Q, ~ 0.1. This 
is particularly the case for the toroidal modes of the model NS13T8. In the same sense, the response of the low radial order 
spheroidal shear modes ;s„ with low I to rotation can be substantial as well, particularly for the model NS13T8. The effects 
of slow rotation on axisymmetric / and p\ modes, however, are not significant, the result of which is consistent with the 
calculation by Saio (1981) for a polytropic star with the index N = 3 and Ti = 5/3, who have however taken into account 
the effects of rotational deformation. 



4 CONCLUSIONS 

We have calculated axisymmetric (m = 0) oscillation modes of neutron star models that have a solid crust and are threaded by 
a dipole magnetic field. We find that the frequencies of the fundamental toroidal torsional modes are not affected significantly 
by the magnetic field as strong as Bs ~ 10 15 G at the surface, and that high radial order torsional and shear modes are 
susceptible to a magnetic field even if the strength at the surface is much less than Bs ~ 10 15 G. Because both spheroidal 
shear and toroidal torsional modes are almost degenerate in oscillation frequency for a given radial order n, the high radial 
order modes with different Is for a given n are easily coupled in the presence of a strong magnetic field. Since the / modes 
(and p modes) are embedded in the sea of high radial order shear modes with various Is that are sensitive to the magnetic 
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field, the / modes suffer mode crossings with the shear modes as Bs varies, and their identity may become ambiguous in the 
presence of a strong magnetic field. 

We find that the g modes in the fluid ocean and the crust/ocean interfacial modes are suppressed in a strong magnetic 
field, the result of which contradicts that obtained by Carroll et al (1986), who showed that the ocean g modes will survive to 
be a magnetic mode, labeled m/g k in their paper, as the magnetic field becomes strong. The reason for the contradiction may 
be attributable to the difference in the geometry of the oscillating stars and of the magnetic field, that is, Carroll et al (1986) 
calculated oscillation modes of cylindrical stars threaded by a uniform magnetic field. The strength of the magnetic field 
necessary to completely suppress the ocean g modes and the crust /ocean interfacial modes depends on neutron star models, 
and if a neutron star has a hot buoyant radiative region in the ocean, the modes can survive a magnetic field of Bs ~ 10 10 G 
or stronger, although their frequency spectrum could be largely modified by the field. 

As a model for burst oscillations observed in many low mass X-ray binary (LMXB) systems (e.g., Strohmayer et al 1997, 
Strohmayer & Bildsten 2004, van der Klis 2004), Heyl (2004) and Lee (2004) proposed that the oscillations are produced 
by low frequency buoyant r modes propagating in the surface fluid ocean of accreting neutron stars in the systems (see also 
Lee & Strohmayer 2005, Heyl 2005). Since no effects of magnetic field are properly taken into account in their analyses, it is 
needed to reexamine the r mode model for burst oscillations, whether the buoyant r modes can survive a strong magnetic 
field of the neutron stars and how their frequency spectrum is modified by the field. 

We have also examined the effects of slow rotation on the axisymmetric oscillation modes, although no effects of rotational 
deformation are included. Since the first order effects of rotation do not appear for axisymmetric modes with m = 0, it is 
the second order effects of rotation that appear first and are due to the Coriolis force when no rotational deformation is 
considered. The second order effects can be important for torsional modes and shear modes with radial order n > 1 in the 
sense that the second order corrections to the frequency due to rotation can be comparable to the frequency itself for Q ~ 0.1. 
But, for SCR 1806-20, for example, the rotation period is estimated to be 7.56s (e.g., Israel et al 2005), implying that the 
underlying object is a very slow rotator in the neutron star standard so that Cl <C 1 and C2O, 2 should be negligible compared 
to wo- The first order correction due to the Coriolis force for non-axisymmetric modes can be found in Lee & Strohmayer 
(1996). 

We find magnetic modes for a strong magnetic field. Here, magnetic modes are regarded as a oscillation mode that exists 
only in the presence of a strong magnetic field, and whose frequency rapidly increases with increasing Bs- We also note that 
the magnetic modes have their oscillation energy dominantly possessed by the magnetic perturbations. The modal properties 
of the magnetic modes found by Carroll et al (1986), labeled m/g k and a k , are not the same as the properties of the magnetic 
modes found in this paper. Note that although the magnetic modes computed by Carroll et al (1986) reside in the fluid ocean, 
the magnetic modes found in this paper have amplitudes both in the fluid ocean and in the solid crust, and a large fraction 
of the oscillation energy resides in the crust, as indicated by Figures 8 and 9. We note that a simple magnetohydrodynamical 
system can support Alfven waves, as well as fast and slow magneto-acoustic waves in fluids (e.g., Sturrock 1994), and we 
may expect almost the same modal structure even for a magnetic solid, in which torsional waves propagate (e.g., Carroll et 
al 1986). At this moment, however, we have no clear classification scheme for the magnetic modes, which makes it difficult to 
obtain a good understanding of the modes. Further studies are definitely necessary for magnatic modes of neutron stars with 
a solid crust, since the magnetic modes could be important observationally for magnetor as an agent triggering instability for 
flares. 

Because of the assumption of non-magnetic core employed in this paper (§2.2), we completely ignore the possible existence 
of magnetosonic modes in the core and the possible coupling between core magnetic modes and crustal shear modes, for 
example. This assumption would be a serious flaw when we are interested in core (or more global) magnetic modes themselves. 
In fact, the crustal magnetic modes found in the present paper could be flawed in the sense that the assumption of non-magnetic 
core excludes from the beginning the possible existence of more global magnetic modes extending form the core to the crust. 
For non-magnetic modes such as crustal shear modes, however, so long as the field strength is less than ~ 10 1B G at the 
surface so that the magnetic pressure is much smaller than the gas pressure (and/or the shear modulus) in the inner crust 
(see, e.g., Fig. 3 of Carroll et al (1986) or Fig. 1 of Piro (2005)), reflection of crustal shear waves at the bottom of the crust is 
effective to establish crustal shear modes well trapped in the crust, and the effects of the coupling with core magnetic modes 
would be minor, for example, on the frequency spectrum of the crustal shear modes. For the field strength much larger than 
10 15 G, the confinement of crustal shear mode amplitudes into the crust would be imperfect and global treatment properly 
including the core will be necessary. As suggested by one of the referees, a differentially rotating magnetor progenitor would 
produce a toroidal magnetic field, for which the complexities possibly caused by mechanical coupling between the crust and 
the core could be avoided. This case may be important and worth careful examination. Recent discussions on the crust-core 
coupling of magnetic modes, leading to more global magnetic modes of neutron stars, are found in Levin (2006), Glampedakis, 
Samuelsson, & Andersson (2006), and Sotani, Kokkotas, & Stergioulas (2006). 

The assumption of non-magnetic core also leads to neglect of another possible role played by perturbed magnetic field 
in the core. Several authors (e.g., Biront et al 1982, Roberts & Soward 1983, Cambell & Papaloizou 1986) have discussed for 
magnetic normal stars that in a deep fluid region where the magnetic pressure is much smaller than the gas pressure, magnetic 
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perturbations become extremely short, decoupled from mechanical perturbations, which leads to a dissipation of the oscillation 
energy. In this paper, however, because of the numerical difficulty in properly treating the abrupt change of wave properties 
between the solid crust and the fluid core, we have employed a simplifying assumption that the fluid core is non-magnetic. 
Although the assumption of non-magnetic core would be reasonably justfied since in the fluid core the magnetic pressure 
Pb is much smaller than the gas pressure p, appropriate estimations of the dissipative effects of extremely short magnetic 
perturbations in the fluid core are necessary to see whether the dissipation is significant enough to damp the excitation of 
the torsional oscillations. We have almost the same difficulty in treating the outer boundary conditions, for which we have 
assumed for simplicity no emission of electromagnetic waves from the surface even if we include the displacement current 
term. Although Carroll et al (1986) (see also McDermott et al 1988) suggested that the damping effects due to emission 
of electromagnetic waves from the surface are negligible except for magnetic modes, labeled m/gk and a*,, fully consistent 
calculations including the damping effects are obviously necessary. 

In the giant flare of SGR 1806-20 observed December 2004, there exist reports on detection of QPOs at 18, 30, and 92.5 
Hz (Israel et al 2005), at 18, 92.5, and 626.5 Hz (Watts & Strohmayer 2006), and at ~90, ~150, 625, and 1835 Hz, and at 
~720 and 2384 Hz but with less significances (Strohmayer & Watts 2006). In the giant flare of SGR 1900+14 observed in 
August 1998, Strohmayer & Watts (2005) have also found detection of QPOs at 53.5, 84, and 155.5 Hz, and at 28 Hz with 
lower significance. Identifying the QPOs of low frequencies Vao <100Hz with the fundamental toroidal torsionl modes with 
various Is is rather secure (McDermott et al 1988, Duncan 1998, Piro 2005), but identification of the middle to high frequency 
QPOs is not straightforward. In fact, QPOs with frequencies ~ 100 to ~ 1000 Hz can be generated by high I' j/io modes, or 
overtone modes ;/i n >i with low I', or low I spheroidal shear modes (S„>i, or core/crust interfacial modes ;S2, depending on 
neutron star models, and there exist no unique solution to identification if the frequency is the only available information. 
For much higher frequency QPOs like that at 2384Hz, low I fundamental modes ;/ must be added to the list of candidates for 
the QPOs, although much more energy will be required to excite the modes to observable amplitudes than that for toroidal 
crustal modes. 

In this paper, we have been only concerned with axisymmetric modes with m = 0. If we extend our calculation to 
non-axisymmetric modes with m / 0, mode coupling effetcs will be much more significant in a strong magnetic field because 
both toroidal torsional and spheroidal shear modes are nearly degenerate in frequency for a given radial order n, which makes 
numerical analysis difficult and tedious. However, together with employing neutron star models constructed with up to date 
equations of state and shear modulus, the extension to non-axisymmetric modes will be inevitable in order to make possible 
serious comparions between theoretical predictions and observations. 



APPENDIX A: OSCILLATION EQUATIONS FOR A MAGNETIZED STAR FOR AXISYMMETRIC 
MODES WITH M = 

In this Appendix, we give the oscillation equations for axisymmetric (m = 0) modes separately for even modes and for odd 
modes. In the followings, for a given matrix F = (Fij), we employ the symbols [F], {F}, and (F) to indicate the matrices 
defined as 

[F] = (F iJ+1 ), {F} = (Fj+ij), (F) = (Fi+i^+i) (Al) 
for i,j = 1, 2, 3, • • •. 

Al Even Modes 

For the solid crust, we employ the dependent variables defined as 

*i = (S,,(r)), z 2 =(ff (r)), zs=(iT l ,(r)), z 4 = a 2 (1 ± (r 3 ^) - [A ] * 2 ) + 2 Ql ^ (r*i) + ^C»b H , (A2) 



P 

z 6 = a 1 (r^-+{l}z 1 )-^-{M 1 }b H , 2 6 = mr^ - ^ [M ] ib T , b H = (i£(r)) , b T = (&£(r)) , (A3) 
where 

lj=2(j-l), l j= 2j, 1$ = 1 + 2(j - 1) for j = l, 2, 3, (A4) 

and the oscillation equations are given by 

dzi 3ri a 2 r . , 1 1 2p s H 
r— — = zi H A 2 2 H 24 C b , (A5) 

dr «3 Q3 Q3 Q3 p 

[M ]r^l = 2 (l- 21) K Zl + ([Mo] \^ [C,]) z 2 ^Kz 4 + ( ±^KC h) b H , (A6) 
dr \ «3 / V I az 1 2Q3 V Q3 P 2 / 
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{M X } = ({M X } - \ {Co}) z 3 - \ib T , (A7) 
r ^T 1 = ( (U - 4 - acu 2 ) VI - — Co(2Qi +Ci) ] zi + ( (V-2ai)[A ] - — Co [QiA + Ci] ) z 2 + ciu>WC «3 



dr \* p 7 V P 



+1/^4 + [Ao] - 4 Ql r^i + ^ (QoAi - Co) b H + Cl u> 2 V (—) * ((QoQi - 1) *i + 2 [Q Ci] z 2 ) , (A8) 
dr p \ c ) 



r ^d7 = y {1} + ~r {Ml} (2<?1 + Zl + y- ciCo ' 2vi - 2ai1 + 2ai < Ao > + {Ml} [QiAo + 2:2 

+ C ia; 2 */l/ {M,}z 3 - {1} z 4 + (V - 3) z 5 + {1} 2air^- + ({Mi} + 1 {Co}) b ff 

ar p V 2 / 

+ Ci lu 2 V (^j (Ao)" 1 (-2 (3Q Qi - 1 + QoCi) zi - [4Q QiA + 8Q Ci] z 2 ) , (A9) 

= -au 2 vVKzi + ciu 2 vV [Mo] z 2 + (-c x d 2 V\ - 2ai(l - ^Ai)) z 3 + (V - 3) z 6 + ([M ] + i [Ci]) ib T 

-ci(i 2 V r (^-) 2 A 1 : 1 (6QiCo+4QiQoA 1 -CiCo)z 3 , (A10) 

where 1 denotes the unit matrix, and 

cu = Lu/y/GM/R' A , v = 2Q/lj, (All) 
and 



ci = .' , v = --—?-, u= — — -, pb = -^t j -, va = \ -^- l A12 

Mr/M dlnr dlnr 8tt V 47r P 

and 

/i 2 4 

qi = — , a2 = Ti — -ai, a3 = ri + -ai, (A13) 

p o o 

and M and 7? are the mass and radius of the star, and G is the gravitational constant. Note that the terms proportional to 
(va/c) 2 come from the displacement current. 

For fluid regions the dependent variables we use are defined as 

and the oscillation equations are given by 

r ^ = (Y- 1 - 3 )y^-r 1 V2 + [Ao]V3 ' (A15) 
r ^iv = (( CltI;2 + rA ^ 1 + {2Ql + 1/1 + (1 " rA ~ u) 2/2 + "pv Co [QlAo + Cl] Va ~ c ^ vCoy * 

+ |Coy 5 + - ci^ 2 (^) 2 ((QoQi - 1) tfl + 2 [Qod] y 3 ) , (A16) 

™ = -5 (rH + m"« + (™ - 5 ™) -3 - ) _1 (A17) 

{M^^W-iw)^!^)" 1 ,,, (A18) 
{Mx} r^Ms = _lP£ { Ml } (2d + Ci) yi - {1} y 2 + f Cl o; 2 l - % {M x } [QiA„ + d]) y 3 - cia> 2 z, {Mx} y 4 



-/ {M,} y 5 - Cl cD 2 (Ao)" 1 (-2 (3QoQi - 1 + Qod) Vl - [4Q QiA + 8Q d] y 3 ) , (A19) 

[M ] = Cl Lu 2 uK yi - Ci uj 2 u [M ] y 3 + Cl cS 2 y 4 - (/ [M ] + \ [d]) y 6 

+Cl ^ 2 ^) 2 AjT 1 (6Q1C0 + 4QiQoAi - C1C0) y 4 , (A20) 
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where 

f = rA-f-+U + 3. (A21) 
1 1 

The non-zero elements of the matrices Co, Cl, K, M , Mi, Q , Qi, A , and Ai that appear in the above equations are 
given by 

(Qo)i,i = Jl+l, {Qo)i+l,i = Jl+2, (Ql)M = JJ+1, (Ql)i,i+1 = Jl+2, (A22) 

(Co)i, i = -(J + 2)J^ 1 , (C )i+i,i = (' + l)J,+ 2 , (Ci)i,i=IJ£i, (Ci) M+ i = -(Z + 3).7£ a , (A23) 

w^ = /ft' (x)< - i+i = -irl' (A24) 

(Ao)i,i = /(! + 1), (Ai)i,i = (i + l)(i + 2), (A25) 

(M )i,i = j-^Jt+i, (Mo) M+ i = (MOi.i = j^fjfti, (Miji+Li^^, (A26) 
where Z = |m| + 2i — 2 for i = 1, 2, 3, ■ • •, and 



(Z + m)(Z - m) 
(2Z-1)(2/ + 1) 



1/2 



(A27) 



A2 Odd Modes 

For the crust, we use the dependent variables defined as 

*i = z 2 = (^(r)), « 3 =(iTj,(r)), 24 = Q 2 (~ (r 3 Zl ) - A z 2 ) + 2ai|; (r«i) + ^ [C ] 6 H , (A28) 

z B =a 1 (r^+z 1 )-^[M 1 ]6« 2 6 = «ir^ - ^ {M } 6* = {bf{r)) , b T = (&?») , (A29) 

where 

li = 2(j - 1) + 1, = 2j, I'j = 2(j - 1) for j = 1, 2, 3, • • • . (A30) 
The oscillation equations are then 

dzi 3ri ,02. 1 1 2p s w 

r— — = zi H A «2 H z 4 Co b , (A31) 

ar as Q3 Q3 Q3 p 

{ M >r£ = 2 (l - 21) { K }Z1 + ( { M 0} - Ig { C l} ) ,2 - JL {*} 24 + W [Co] - II) 6-, (A32) 

[Mr] = ([Mr] - \ [Co]) 23 - \ib T , (A33) 

= ~ 4 ~ Cl<Ii2 ) V1 ~ ^Y Co (2Ql + 21 + ( (v " 2ai) A " " "T^ (<5lAo + 22 + Pol 23 

+V24 + A z 5 - 4 Ql r^i + [Q Ai - Co] 6 H + ciu>V ( — Y ((QoQi - 1) zi + 2Q C lZ2 ) , (A34) 
ar p \ c J 

= ( V1 + V" Ml (2Ql + 21 + ( K - C1L ° 2V1 - 2ai ( X " A( ») + ^T" Ml ( QlAo + Z2 

+ Ci lu 2 vV [Mi] 23 - Zi + (V - 3) 25 + 2air^- + ([Mr] + i [C ]) b H 

+cia) 2 V (^-) 2 A _1 (-2 (3Q Qi - 1 + Q0C1) zi - (4Q QiA + 8Q0C1) z 2 ) , (A35) 

r ^££ = -ciwW {.ftT} zi + ciu> 2 z/V {Mo} 22 + (-ciwVl - 2ail + ai (Ai)) 23 + (V - 3) z 6 
ar ' v ' 

({M }+ i{Cr}) ib T -ciciV (^-) 2 (Ai)" 1 [6QiCo + 4QiQoAi -CiCo]23. (A36) 
For fluid regions the dependent variables we use are 
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( P ' l i ^ A 4 P B uH ^PB ., T , a ovN 

Vi = zi, y 2 =1-2^1, Vi = z 2, Va = zz, y 5 = —b , y 6 = —ib, (A37) 
and the oscillation equations are then given by 

r ^ = (f- 1 - 3 )y^-f 1 y2+AoV3 ' (A38) 

r^l = ^( Cl cu 2 + tA) 1 + 5* C„ (2Q 1 + d)j Vl + (1 - rA -U)y 2 + ^C Q (QiA + d) y 3 - Cl u) 2 zv [C„] y 4 

+^ Po] y 5 + ~ ClC ° 2 (t) 2 ((0(,Ql " 1} Vl + 2QoCiy ^ ' ( A39 ^ 

W ^ = -K^- 4 ){^ + ^ W + (< M °> - 5 ™) »» - 2" (>) " y5 ' (A40) 

^t<^->^-i{wy ly - (A4i) 

dy 5 _ 4p s R/f . i ^ i ^___ 2l 8p s 



[Mi] r-^- = -^Mi (2Qi + Ci) yi - y 2 + ^ciO'l - ^Mi(QiA + Ci)J y 3 - cxtfv [Mi] y 4 - / [Mi] y 5 

-cia. 2 2 Aq- 1 (-2 (3Q Qi - 1 + QoCi) yi - (4Q QiA,, + 8Q Ci) y 3 ) , (A42) 

{Mo} = cio/V {#} Wl - cio; 2 ^ {M } y 3 + ci^ 2 y 4 - (/ {M„} + 1 {Ci}) y 6 

+cio; 2 (2±) 2 (Ai)- 1 [6QiC + 4QiQ Ai - C1C0] y 4 . (A43) 
The non-zero elements of the matrices Co, Ci, if, Mo, Mi, Qo, Qi, Ao, and Ai are given by 

(Qo)i,i = JT , (Qo)i,i+l = Jl+l, (Ql)i,i = JF i (Ql)i+l,i = Jl+1, (A44) 

(C )i,i = (l-l)Jr, (Co)i,i+i = -(l + 2)JlZ 1 , (C 1 )i, i = -(l + l)JT, (Ci)i, i+ i = {iW+i, (A45) 

jm Jin 

(K)i,i = -±-, (JOm-h = jji-, (A46) 
(A )i,i = /(/ + l), (Ai) M = (/ - 1)1, (A47) 

(M ) M = ^jr, (Mo) i+ i,i = j^yJffi, (Mihi = ^pJT, (Mi), + i,^Jm, (A48) 

where Z = |m| + 2i — 1 for i = 1, 2, 3, • • •. 

Note that for axisymmetric modes with m = 0, the components H tl and 6^ for even modes and iT t i and 6^? for odd 
modes are missing. As is evident from the equations given above, if Q — 0, the oscillation equations for axisymmetric modes 
are separated into those for spheroidal and toroidal modes, respectively. In this case, we can obtain spheroidal modes and 
toroidal modes separately, and mode couplings between these two are brought about as a result of the effects of rotation. 



APPENDIX B: BOUNDARY CONDITONS AND JUMP CONDITIONS 

The surface boundary conditions we use are the same as those given in Lee (2005), and they are at the stellar surface 

Sp/p = 0, b s + L + b H = 0, (Bl) 

and 

ib T = 0, (B2) 
where (L + )ij — 5^(1'^ + 1), 8p is the Lagrangian perturbation of the pressure, and 

b s = (bf, (r)) = -A 1 Kzi - 2AiM Z2 - 2mz 3 . (B3) 

3 

Assuming the central region is a non-magnetic fluid core, the central boundary conditions we impose are the regularity 
conditions for the variables ry 1 and ry 2 at the center (see, e.g., Lee & Saio 1986). 

For the jump conditions at the interfaces between the solid crust and fluid regions, we assume that the dipole magnetic 
fiels is continuous at the interfaces, and that there appears no singularity in B' at the interfaces. Then, with the continuity 
of the radial component of the displacement vector given by 
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zi = Vi, 

the conditions 

[(£ x B ) ■ e e]t = 0, 
lead to 

= y 3 , 
and 

23 = y 4 , 



[(I x B) ■ e^]+ = 0, 



(B4) 



(B5) 



(B6) 



(B7) 



where [F]^_ = lim e ^o [F(r + e) — F(r — e)]. Note that the conditions (B5) are obtained by integrating the right hand side 
of equation (5) along a closed, infinitisimally small, curve that crosses the interface. The conditions (B4), (B6), and (B7) 
guarantee the continuity of the radial component of B' at the interfaces because of equation (-B3). Finally, we require the 
continuity of the three components of traction given by 

\5T^ r + + = 0, 



and 



0, 



(B8) 
(B9) 

(BIO) 



[Sr^f, + St^] * = 0, 

where St? and 5t? are the ij components of the perturbed traction associated with bulk modulus and magnetic field, 
respectively, and they are given by 



St": 



dY? 



Sr r<t> — ai 



E{ Q2 [^Jr^O-'C + l)".] +2a 1 |(rS I )}y« m , 
i 



+ r- 



dT v 1 



dY™ 



dr sin 



i,v 



1 917; 
in d<j> 



dr 86 



5r rr 
and 



- — (B r 5B r 

47T 



B e 8Bg) , 



— {BgSB r + B r SB e ) , 

47T 



5r, 



47T 



-B r SBj, 



Note that 

[Bi^-Ji = [Si (B; + 1 ■ VB 3 )] + = [BiB' 3 ] 



(Bll) 
(B12) 

(B13) 

(B14) 
(B15) 

(B16) 

(BIT) 



since =0 and [B]^ = at the interface. 

In this paper, all the jump conditions discussed above are applied at the interface between the solid crust and the fluid 
ocean. However, since we have assumed the fluid core extending below the solid crust is non-magnetic, we cannot use all the 
jump conditions at the interface between the solid crust and the fluid core. Here, we apply the jump conditions (B4), (B8), 
(B9), and (B10). In the limit of non- magnetized stars, these jump conditions reduce to free slippery conditions employed by 
McDermott et al (1988), and Lee & Strohmayer (1996). Because of the jump conditions we use, the horizontal and toroidal 
components of the displacement vector and the perturbed magnetic field are not necessarily continuous at the interface. 
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Figure 1. u> versus log(Bs) for axisymmetric (m = 0) toroidal torsional modes ;/t n of the neutron star models NS05T7 (top panel) 
and NS13T8 (bottom panel), where the fundamental torsional modes ; /t n= o for V = 2 to 5, and the overtone modes ;/_ 2 t n for n = 1 to 

4 for NT05T7 and for n = 1 and 2 for NS13T8 arc shown. The frequency Q shown is the local one, and i/oo = w (l - 2GM/Rc 2 ) 1/2 /2tt 
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Figure 2. Eigenfunctions iTy and ibj of the axisymmetric (m = 0) fundamental toroidal torsional mode 2*0 f° r Eg = 10 12 G in panels 
(a) and (b) and for Bg = 10 15 G in panels (c) and (d), where the solid lines stand for iT;/_ 2 and ibJL-^, and the dashed lines for iTy—^ 
and ibT_ 3 , respectively. The amplitude normalization is given by iT ( /_ 2 = 1 at the surface. 



Axisymmetric oscillations of magnetic neutron stars 17 



CM 



CO 

DC 



0.1 



0.01 




1000 



N 



1 00 > 



12 13 14 

Log (B ) 



CO 



0.1 



0.01 




1000 



N 



1 00 



Log (B ) 

a icr s' 



Figure 3. ui versus log(B$) for axisymmetric (m = 0) spheroidal shear modes is n and core/crust interfacial modes ; «2 of the neutron star 
models NS05T7 (top panel) and NS13T8 (bottom panel), where the interfacial modes ;«2 have large amplitudes at the interface between 
the solid crust and the fluid core. For the model NS05T7, magnetic modes, expediently labelled mj,, are also plotted. The frequency u> 
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and for B$ = 10 15 G in panels (c) and (d), where the solid lines stand for Hi =2 and bff_ 1 , and the dashed lines for i?; = 4 and bf}_ 3 , 
respectively. The amplitude normalization is given by 5; = o = 1 at the surface. 
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Figure 5. w versus log(Bs) for axisymmetric (m, = 0) fundamental 1—2/ mode of the neutron star model NS05T7. The i—^f mode 
suffers mode crossings with high radial order shear modes ;s„ as Bg increases. The frequency u) shown is the local one, and Uac = 

u> (l — 2GM/Rc 2 ^j ^ /2-7T is also dsiplaycd on the right axis of the panel. 
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Figure 6. w versus log(B 3 ) for the axisymmetric (m = 0) ;=29n=l mode and the crust/ocean interfacial mode ;=2*l whose amplitudes 
peak at the crust/ocean interface for the neutron star model NS05T7. 
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Figure 7. ui versus log(B 3 ) for the fundamental torsional mode 2*0 (left panel) and the first overtone torsional mode 2*1 (right panel) 
of the model NS13T8. The panels are magnifications of the corresponding parts from Figure 1. 
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Figure 8. Eigenfunctions iT v and ibj of the m = toroidal magnetic mode of Q = 0.01557 at B s = 10 14 ' 48 G for the model NS13T8, 
where the solid lines stand for iT;/ =2 and ibf =1 , and the dashed lines for iT ( / =4 and ibJL 3 , respectively. The amplitude normalization is 
given by iT ; / =2 = 1 at the surface. 
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Figure 9. Eigenfunctions Hi and b^j of the m = spheroidal magnetic mode labeled mi at Bg = 10 14,45 G for the model NS05T7, 
where the solid lines stand for i?; = 2 and b^_ 1 , the dashed lines for i?;=4 and bfi_ 3 , and the dotted lines for Hi=§ and b^f_ 5 , respectively. 
The amplitude normalization is given by Si = q = 1 at the surface. 



